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We use the marginal stability condition for galactic disks and the stellar ve- 
locity dispersion data published by different authors to place upper limits on 
the disk local surface density at two radial scalelengths R = 2h. Extrapo- 
lating these estimates, we constrain the total mass of the disks and compare 
these estimates to those based on the photometry and color of stellar popu- 

*5 lations. 
i The comparison reveals that the stellar disks of most of spiral galaxies 

2 in our sample cannot be substantially overheated and are therefore unlikely 

to have experienced a significant merging event in their history. The same 

l^j conclusion applies to some, but not all of the SO galaxies we consider. How- 

ever, a substantial part of the early type galaxies do show the stellar velocity 
j> dispersion well in excess of the gravitational stability threshold suggesting a 

major merger event in the past. We find dynamically overheated disks among 

,_i both seemingly isolated galaxies and those forming pairs. The ratio of the 

marginal stability disk mass estimate to the total galaxy mass within four 
radial scalelengths remains within a range of 0.4 — 0.8. We see no evidence 
for a noticeable running of this ratio with either the morphological type or 

t-h color index. 

J^j 1. Introduction 

The majority of stars in spiral and SO galaxies are in general concentrated in 
old stellar disks with characteristic half-width of several hundreds parsecs. 
A formation and dynamical evolution of disks are the matter of hot debates, 
so the analysis of disk kinematical characteristics plays a key role (see e.g. 
Zasov, Silchenko (2010) and references therein). By kinematical characteris- 
tics we mean a rotation velocity and velocity dispersion of old stars, making 
up the bulk of stellar disk mass, as functions of radial coordinate in a disk 
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plane. Typical velocity dispersion of old stellar disks equals to several dozens 
km/s. In general, it either can reflect the velocity of turbulent motion of gas 
which gave the birth to the first stellar generation, or can be the result of 
the dynamical heating of disk during its evolution. 

Several possible mechanisms of equilibrium disk heating are known, the 
most effective are a scattering of stars on massive clouds or stochastic spirals, 
gravitational interaction with dark halo subunits, major and minor mergers 
(see Jenkins, Binney (1990); Binney, Tremaine (2008); Shapiro et al. (2003)). 
As far as the velocity dispersion of stars increases the efficiency of all mech- 
anisms decreases. For initially dynamically cold disk, another mechanism 
capable of effectively heating a disk is the development of gravitational in- 
stability against random perturbations. It leads to the increase of stellar 
velocity dispersion both in the plane of a disk, and along the vertical axis 
(through the developing of bending instability). At the early stage of evo- 
lution during the disk mass growth, stellar velocity dispersion might remain 
fitted to the level of its marginal stability. If a total disk mass is constant 
and there exists a gas disk, forming stars with low velocity dispersion in its 
plane, a self-regulating process, supporting stellar disk velocity dispersion 
in the quasi-stationary station, is also possible (Betrin, Lin, 1987). In any 
case, even if gas and young stars are absent, the minimal velocity dispersion 
of old stars at given R will be constrained by disk local marginal stability 
condition. Numerical simulations of models of initially weakly unstable 3D 
disks show the fast passage of a disk into the marginally stable state, after 
which the increasing of velocity dispersion practically ceases in the absence 
of the heating mechanisms listed above (see e.g. Khoperskov et al. (2003)). 

For a thin disk with parameters slowly varying with R the local critical 
value of radial velocity dispersion of stars is defined by Toomre criterion: 
c r = ct, where ct = 3.36Gcr / k, a is the disk surface density at given R, k is 
epicyclic frequency determined by angular velocity Q and its radial derivative: 

k = IVl ■ yJ(l + (R/2Sl) ■ (dSl/dR)). (1) 

A finite disk thickness makes disk more stable, while non-radial perturbations 
have the opposite effect (they decrease the stability threshold). The stability 
condition taking into account both of these effects can not be expressed in 
analytic form. If to express the critical radial velocity dispersion as (c r ) cr = 
Q-ct, then, as numerical simulations of 2D and 3D disks show, the parameter 
Q for a wide range of R lies in the interval 1.2 — 2.5 (see e.g. Bottema (1993), 
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Khoperskov et al. (2003)). Note, that the analytically obtained local stability 
condition of 2D disk taking into account non-radial perturbations gives for a 
flat rotation curve Q c ~ 3 (Polyachenko et al. (1997)). 

The assumption of marginally stable disks makes it possible to put con- 
straints on its surface density: 

a{R) ~ 3.36GQ(RY (2) 

and hence on a total disk mass therefore. It enables to obtain relative masses 
of disk and spheroidal components (bulge and dark halo) (Bottema (1993), 
Bottema (1997), Tyurina et al. (2001), Zasov et al. (2006)) and to estimate 
local thickness of stellar disk (see f.e. Khoperskov et al., 2010), which allows 
in turn to find a volume density of gas layer in disk plane if it exists in a 
galaxy (Zasov, Abramova (2006), Kasparova, Zasov (2008), Zasov, Abramova 

(2008) ). 

However the question of how far are the real disks from marginally stable 
state still remains open. There are two indirect evidences of closeness of stel- 
lar velocity dispersion to (c r ) cr for significant fraction of disk galaxies. First, 
disk mass-to-light ratios, which follow from the stability condition are in good 
agreement with the estimates based on stellar population evolution photo- 
metrical models (see e.g. Bottema (1997), Zasov et al. (2004)). Second, the 
assumption of disk marginal stability enables to explain a positive correlation 
of relative disks thicknesses and their relative masses within the optical bor- 
ders (Zasov et al. (2002)). Note, however, that both numerical simulations 
(Sotnikova, Rodionov (2006), Khoperskov et al. (2010)) and the analysis of 
observations of edge-on galaxies (Mosenkov et al. (2010), Bizyaev, Mitronova 

(2009) , Bizyaev (2010)) show that this correlation is rather sparse especially 
for galaxies with small bulges. Moreover, numerical dynamical models by 
far are not always compatible with the assumption of disk marginal stability 
(see, e.g., Zasov et al. (2008)). 

A difficulty of estimation of the ratio of observed velocity dispersion c Q b s 
to its critical value (c r ) cr lies first of all in the difficulty of measurement of 
c b s at sufficiently large galactocentric distances and in addition is due to a 
number of factors hard to be taken into account even in the axissymmetrical 
galaxy model. Among the latter factors are the variation of critical Toomre 
parameter with radius, the ratio of radial to vertical velocity dispersion and 
the account of asymmetric drift in the estimation of disk angular velocity 
especially for the early type galaxies. 
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The main goal of this work is to find out how far the observed values 
of velocity dispersion at given R are from the expected ones for marginally 
stable disks for different type galaxies. For this purpose, we analyze the 
velocity dispersion data, available in the literature, for a chosen galactocentric 
distance R = 2h, where h is the radial disk scalelenght. At this radius a 
disk contribution to the rotation curve is nearly maximal. At higher R, 
measurements of velocity dispersion are less reliable and also the influence of 
external factors dynamically heating a disk may be more important. In turn, 
at lower R it is often hard to separate the contributions of bulge and disk 
stars to the measured values of c & s , and to take into account the presence 
of a bar if it occurs in a galaxy. It is worth noting that, as numerical 3D 
N-body simulations show for different parameters of components (disk, bulge 
and halo), the spread of values of Q(R) is minimal at R ~ 2h. Following 
Khoperskov et al. (2003) here we assume Q(2h) ~ 1.5. 



2. The sample and the method we use 



To estimate the upper limits of disk masses corresponding to the marginal 
stability condition we have chosen objects for which sufficiently extensive 
radial distributions of stellar velocity dispersion and rotation curves were 
available in the literature (121 objects). Using these data we estimated the 
threshold value of surface density corresponding to the gravitational stability 
at R = 2h from the equation (2). Radial velocity dispersion of stars c r , if it 
was not given in the cited literature, was determined from c f, s (R) measured 
along the major axis of a galaxy: 



^obs 



2 • 2 • , 2 2 • 

= Cj. sin % + c z cos % 



(3) 



which leads to: 



Cobs 



(c^/cr) 2 sin 2 i + (c z /c r ) 2 cos 2 



-1/2 



(4) 



where and c z are azimuthal and vertical components of velocity dispersions 
correspondingly, and i is the inclination of a disk. Radial and azimuthal 
velocity dispersions are interconnected by Lindblad formula c^/c,. = k/2Q 
following from the epicyclic approximation. Direct measurements of c z /c r 
ratio in galaxies show that in majority of cases it lies in the interval 0.4 — 0.7 
(Shapiro et al. (2003)). In current work it was taken to be 1/2. 
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In most cases rotation curve V(R) has a plateau at the distance of two disk 
radial scalelenghts, so the k/VL ratio can be taken as \/2. For the increasing or 
decreasing rotation curves we applied the corresponding corrections. If only 
the stellar rotation curve K(-R) was available for a galaxy, the approximate 
asymmetric drift correction was introduced following Neistein et al. (1999). 
Note, that for c r /V c > 0.5 this correction is inevitably very crude. 

A list of galaxies with available measurements of stellar disk velocity 
dispersion and their observed and calculated properties is given in Table 
1. For galaxies with line-of-sight velocity exceeding 800 km/s the adopted 
distances correspond to the Hubble constant H = 75 km/s/Mpc, (except 
Virgo cluster galaxies for which we took D — 17 Mpc). For several nearby 
galaxies with lower line-of-sight velocity (NGC 3198, NGC 6503, IC750) we 
used the distance modulus given in Hyperleda database. The disk inclination 
i, color indices corrected for galactic extinction, inclination and redshift, and 
line-of-sight velocities were also taken from this database. For the cases when 
the inclinations taken from Hyperleda differed significantly from those given 
in the original sources referred in the Table, the preference was given to the 
later ones. 

After the estimation of surface density at R = 2h, it is easy to find a total 
disk mass: 

M d = 2nh 2 a(2h) e 2 . (5) 

In this work a disk photometric radial scalelenght and the scalelenght of 
surface density are assumed to be equal. In some cases this can be quite rough 
approximation, but the introduced uncertainty is not so high in comparison 
with other error sources: a change of h by ± 30% for a fixed adopted value 
a(2h) leads to the change of M d by +7%, -15%. 

To compare the obtained disk masses with the photometric estimates 
we calculated disk mass-to-light ratios in B band (M/L B )d- For Sb and 
later type galaxies where the bulge contribution to the luminosity is low, 
we neglected the difference between the total and disk luminosity. For Sa — 
SO galaxies we have taken into account the contribution of bulge. To do 
this we took a disk luminosity estimates or disk de-projected surface bright- 
ness at R — 2h from de-projected radial surface brightness profiles avail- 
able in the literature. In the absence of such data a disk luminosity was 
calculated from the total luminosity using bulge-to-disk luminosity ratio 
B/D, average for given morphological type (taken from Graham (2001)): 
L d = L tot /((B/D) + 1). Disk color index was obtained in this case from the 
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total color index: 

(B - V) d = 2.5 • lg (l0(° A <B-v) tot ) _ B / T . 10 (,o.HB-v) bulgc ^ _ (6) 

Bulge color index (B — V\ u i g( , of early type galaxy was taken to be 0.8, if 
the total color index of galaxy did not exceed 0.8. Otherwise (for very red 
galaxies), a color index of a disk was accepted to be equal to the total one. 

3. Results 

Surface density and disk mass estimates together with other data and the 
references to the data sources are listed in Table 1. It contains: 

(1) Name. 

(2) Morphological type. 

(3) Radial disk scalelength h used in the paper. 

(4) Ratio of radial velocity dispersion to circular velocity at R = 2h. 

(5) Disk surface density corresponding to the marginal stability condition at 
R = 2h. 

(6) Mass-to-light ratio of a disk, marginally stable at R = 2h (5-band). 

(7) Ratio of mass of the disk, marginally stable at R = 2h, to the total mass 
M t within R = 4h, where M t = 4V^h/G and V c is circular velocity. Cases of 
Md/M t > 1 indicate a significant dynamical overheating of a disk. 

(8) Velocity dispersion and rotation velocity data sources. 

Diagrams u {M/Lb)ci — {B — V)o" for stellar disks are shown in Fig. 1, 
where disks masses M d are obtained from the gravitational stability condition 
at R = 2h^ For the Galaxy the mean color index for Sbc type was used 
(Buta et a!., 1994). A straight line reproduces model relation obtained by 
Bell, de Jong (2001) for old stellar systems with different present-day star 
formation rate (SFR) using modified (bottom-light) Salpeter initial mass 
function (IMF). It is worth noting that the light extinction inside the galaxies 
that may play important role for disks with high content of dust decreases the 
luminosity of a galaxy and simultaneously increases its color index shifting 
the points practically along the theoretical correlation. 

In Fig. la the entire sample of galaxies with known color indices is shown; 
pair members are marked by asterisks. Here and in the other figures the er- 
ror bars correspond to the errors of velocity dispersion used to calculate the 

1 Color indices are available not for all considered galaxies, so the number of points in 
the diagram is not equal to the number of galaxies listed in Table 1. 
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disk local surface density. The uncertainty of parameter Q c is not included 
in the error bars (according to the numerical simulations it reaches at least 
20%). The same diagram as in Fig. la, but after the exclusion of SO- a — SO 
galaxies, is shown in Fig. lb. Though the scatter of points in these dia- 
grams is large (~ 0.3 dex), it is compatible with the errors of individual 
mass estimates. Therefore one can conclude that there is a general agree- 
ment between (M/Lb)<i estimates based on marginal stability condition and 
those based on the stellar population evolution modeling. It is remarkable 
that most of the galaxies, which significantly deviate from model relation 
have a red color (B — V) > 0.7. At least half of these galaxies are above the 
straight line, that is most of them have disks with significant dynamical over- 
heating. Note, that there are both paired and isolated systems among these 
galaxies. Pair members occur also among the late type galaxies in Fig. lb 
(11 objects), but the presented statistical data are not enough to claim their 
systematical difference from the isolated systems. Anyway, pair, group and 
cluster membership is not always connected with the stellar velocity disper- 
sion exceeding the level needed for marginal gravitational stability. It seems 
that only a strong gravitation perturbation can be a cause of disk dynami- 
cal overheating. This conclusion also follows, for example, from the higher 
mean thickness of stellar disks in interacting edge-on galaxies (Reshetnikov, 
Combes (1997)). 

The position of our Galaxy in the diagrams is shown by the open circle. 
For our Galaxy {M/Lb)ci ratio was calculated using the local values of sur- 
face density and brightness measured for solar vicinity and accepting h = 3 
kpc. The surface brightness data in 5-band was taken from Portinari et al. 
(2005).. Note, that the value of surface density in the Solar neighborhood 
obtained on the basis of marginal stability condition is in a good agreement 
with the direct estimates (Korchagin et al. (2003), Holmberg, Flynn (2004), 
Kuijken, Gilmore (1991)). This gives evidences of marginal stability of disk 
of Galaxy at least at the galactocentric distance of a few radial scalelenghts. 
This conclusion conforms with the idea of the absence of major merger events 
in the history of our Galaxy (see e.g. Wyse, 2009)). 

In Fig. 2a the ratio of radial velocity dispersion to circular velocity c r /V c 
is correlated with morphological type of galaxies. Type t is coded as de 
Vaucouleurs numerical type: spiral galaxies correspond to the interval t=l 
(Sa) — 7 (Sd), higher values of t represent irregular galaxies, and t = — 1...— 3 
are related to SO galaxies. The correlation between c r /V c and type is rather 
weak, moreover, if to exclude SO galaxies, it almost disappears. The same 
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conclusion may be done for the correlation between c r and V c at the fixed 
radius R = 2h (Fig. 2b). It means that it is impossible to get a reliable 
estimation of disk velocity dispersion by the indirect way, from the type of a 
galaxy or its rotation velocity, as it was for example suggested by Bottema 
(1993). 

The conclusion that most of the galaxies with dynamically overheated 
disks are among the red galaxies (Fig. la) is confirmed by Fig. 3, where 
color index is compared with the relative disk mass Md/Mt within the optical 
radius R = Ah. As above, here Md corresponds to disks, which are assumed 
to be marginally stable at least at R = 2h, and total mass of a galaxy was 
calculated as M t = Ah ■ V?/G. 

As far as Md we found is in general case the upper limit of mass of a 
disk, the values Md/Mt > 1 obtained for some galaxies (see Fig. 3) have no 
physical sense and indicate unambiguously the significant overestimation of 
disk masses, or, in other words, the disk overheating. Almost all of these 
objects are lenticular galaxies with a high color indexes (B — V)q > 0.7. 

Lenticular galaxies differ from other galaxies of the sample by their dis- 
tributions of both Md/M t ratios (Fig. 3) and "marginal" surface densities 
at R = 2h (Fig. 4), being shifted to higher values of Md/M t and densities. 
These differences are statistically significant: according to the Wilcoxon cri- 
terion, the level of significance a = 0.99. From the distribution of Md/Mt 
(Fig. 3) it follows that for most of galaxies disk to total mass (inside R = Ah) 
ratios lay in the interval 0.4 — 0.8, independently of color and morphological 
type (if to exclude galaxies with apparently overheated disks, for which these 
ratios remain uncertain). 

If to consider the disky galaxies as systems with marginally stable disks, 
one can plot the relation Md — V c (baryonic Tully-Fisher relation), connecting 
the most important parameter of stellar disk (a mass) with the parameter de- 
termined mostly by massive halo (rotation velocity). This diagram is shown 
in Fig. 5 (galaxies with Md/M t > 1 are excluded). The relation obtained by 
McGaugh (2005) for a large sample of galaxies, disk masses of which were 
estimated from their luminosities and colors, is also shown by straight line 
for comparison. From this figure it follows that disk masses estimated from 
the marginal stability condition show the same relation as masses obtained 

2 Note that this simple formula for a total mass, is, strictly speaking, correct for spheri- 
cally symmetric distribution of density and may underestimate Md/M t by 20-30% Md/M t 
for a disk galaxy which does not possesses a massive bulge or halo. 
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from photometry and evolution models, which confirms the absence of sys- 
tematical difference between these evaluations. 

4. Conclusions 

As it was shown above, if to assume the marginal stability of inner parts of 
disks of spiral galaxies, one can estimate local values of disk surface densities 
and as the next step — the absolute and relative values of their total masses 
or, strictly speaking, the upper limits of masses for stable equilibrium disks. 
These estimates in most cases appear to be in a good agreement with the 
values obtained from photometric data on the basis of stellar population 
evolution models (Fig. lb, Fig. 5). It allows to admit that the disks of most 
spiral galaxies are close to the marginally stable state. 

Evidently, mass and density estimates presented in this paper for any 
given galaxy remain rather rough, mainly due to large errors of radial ve- 
locity dispersion c r , and, to a lesser degree, of radial scalelength h values. 
For some of the galaxies the resulting errors exceed a factor of 2 (see Fig. 
1). Thus to obtain more reliable estimates for individual galaxies one have 
to construct dynamical models consisted of several components which repro- 
duce the radial profiles of observed rotation velocity and velocity dispersion 
c bs of disk stars. Still, the estimation of disk masses for galaxies of different 
types presented in this paper allows to conclude that there is no significant 
systemic dynamical overheating of disks in most of galaxies, except lentic- 
ulars, where the velocity dispersion often exceeds the threshold needed for 
marginal stability. It follows that many spiral galaxies have not suffered a 
strong gravitation perturbations or mergings after the formation of bulk of 
their stars. This result creates a certain problem in the generally accepted 
theory of galaxies formation in a gravitational well of dark halo (see e.g. the 
discussion in Zasov, Silchenko, 2010). The conclusion of marginal stability 
of disks, however, may not concern their outermost regions, which require an 
individual analysis and often may be dynamically overheated. 

It is essential, that the excessively high velocity dispersion at R — 2h 
usually takes place in galaxies which, if to judge by color indexes, have a 
weak star formation (color index (B — V) > 0.7). Most of them belong to 
SO /a — SO types. For these galaxies almost total absence of young stars and 
gas parallel with high stellar velocity dispersion apparently resulted from at 
least one merger with massive satellite in their history. This event could both 
increase the velocity dispersion of stars and decrease the content of cold gas 
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evidently as the result of intensive star formation burst triggered by merging. 
Gas partially may transform into stars, partially could be ejected from the 
disk due to the activity of a large number of massive stars. 

Nevertheless, as it follows from Figs la and 3, some "red" galaxies still 
appear to keep their disks in the state close to marginal stability despite a 
high color index. It supports the idea of the existence of different mecha- 
nisms responsible for the transformation of spiral galaxies into lenticulars. 

This work was supported by grants RFBR 07-02-00792, 09-02-97021. 
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Table 1: Main results 
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Figure 1: Diagram "(M/L^)d — {B — V)q \ plotted for (a) the entire sample of 
galaxies and (b) for the galaxies of types later than SO /a. Milky Way is marked 
by open circle, pair members are shown by asterisks. Straight line corresponds to 
model relation obtained for old stellar systems with different history of star forma- 
tion by Bell, de Jong (2001). Here and in other figures the error bars correspond 
to the errors of velocity dispersion data used to estimate the disk masses. 
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Figure 2: (a) Diagram U c r /V c — morphological type", (b) Diagram U V C — c r ". 
Galaxies with types later than Sa are marked by open circles; a straight line 
corresponds to the mean value of c r /V c . 



22 



1.5 



O 



1,0 • 



O DO V 

o o _V * 



o 



o 



o 



0,5- 



o 



o 



o 



o 



Q> 



o 

0°° 

o 



o ~ o 



CP 



o 



o 



o 



0,0 • 



— I ' 1 1 1 1 1 — 

0,6 0.7 0,8 0,9 



— I 

1,0 



0,3 0,4 



0,5 



(B-V)„ 



Figure 3: Diagram "M d /M t — (B - V) " ■ A horizontal line marks M d /M t = 1. 
Points above this line correspond to the galaxies with certainly overheated disks. 
Open circles represent the galaxies with types later than Sa. Typical errors of 
M d /Mt caused by the uncertainty of stellar velocity dispersion estimates are about 
40%. 
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Figure 4: Diagram u a(2h) — morphological type". Pair members are marked by 
asterisks. Straight horizontal line separates galaxies with a(2h) = 500M sun /pc 2 
level, which definitely evidences the overestimation of their mass. 
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Figure 5: Baryonic Tully-Fisher relation between the disk mass and rotation ve- 
locity. Straight line corresponds to the relation presented by McGaugh (2005) for 
spiral galaxies where disks masses were found from their luminosity and color . 
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